o 



Mon. Not. R. Astron. Soc. 000, [THS] (2010) Printed 26 July 2010 (MN IM^X style file v2.2) 



Accretion disk in the eclipsing binary AU Mon 



G. Djurasevic 1 ' 2 !, O. Latkovic 1 , I. Vince 1 and A. Cseki 1 

1 Astronomical Observatory, Volgina 7, 11060 Belgrade, Serbia 

2 Isaac Newton Institute of Chile, Yugoslavia Branch 



Accepted 2010 July 6. Received 2010 July 4; in original form 2010 April 12 



6 



ABSTRACT 

We analyze the CoRoT and V-passband ground-based light curves of the interacting 
close binary AU Mon, assuming that there is a geometrically and optically thick ac- 
cretion disk around the hotter and more massive star, as inferred from photometric 
and spectroscopic characteristics of the binary. Our model fits the observations very 
well and provides estimates for the orbital elements and physical parameters of the 
components and of the accretion disk. 



Key words: 

cretion disks. 



stars: binaries: eclipsing - stars: individual: AU Mon - accretion - ac- 



> 

(N 

oo 
m 

o 
o 



% 



1 INTRODUCTION 

AU Mon (HD 50846, HIP 33237) is an eclipsing, double- 
lined spectroscopic binary. The system consists of a Be-type 
primary (which will be referred to as the gainer in the rest 
of the text), and an evolved G-type secondary (which we 
will refer to as the donor), that has likely filled its Roche 
lobe and is losing mass to its companion. Spectroscopic evi- 
dence indicates the presence of several layers of circumstellar 
matter in the system, a nd allows for the existen ce of a per- 
manent accretion disk (ISahade fe Ferrer 1 1 19821) . According 
to these characteristics, AU Mon belongs to the class of hot, 
massive Algols. 

The orbit of the system is circular, with an orbital pe- 
riod of about 11 days. There is an additional periodicity in 
the total light of the system, with a period of about 417 days 
and an ampli tude of about 0.2 5 mag. This variabi lity was 
discovered b y lLorenzi I (I1980al ) and interpreted by IPeters I 
l|199ll . [199J) as resulting from periodic fluctuations in the 
mass transfer rate, which may be caused by oscillations of 
the donor star about its critical Roche surface. 

A thorough overview of previous knowl edge about the 
system can be found in the recent paper by iDesmet et al. I 
(J2010r) , who performed a very detailed analysis of the system 
based on C0R0T photometry (the same data this study is 
based on) , a collection of previously published ground-based 
photometry, and on high-resolution spectroscopy. Based on 
the analysis of ground-based data, they concluded that the 
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long-term variation originates from changes in the trans- 
parency of circumbinary material. They have also found pe- 
riodicities shorter than the orbital period in C0R0T data. 
Two reliable frequencies of 10.4 and 8.3 d _1 were tentatively 
assigned to the pulsation of the gainer, and the power excess 
in low-frequency region to the solar-like oscillations of the 
donor. 



However, the attempt bv lDesmet et al. I (|2010l ) to model 
AU Mon as a semidetached system wasn't entirely suc- 
cessful. Most notably, the semidetached model could not 
repro duce the Ros siter-McLaughlin effect (IRossiter I Il924l ; 
IMcLaughlin 1 1 19241 ) . and required an unusually large grav- 
ity darkening exponent to fit the light-curves. The authors 
therefore suggested that the system should be analyzed us- 
ing a model which includes the gas stream and the accretion 
disk around the primary component. 

There are several spectroscopic indications for the ex- 
ist ence of an accretion d isk in AU Mon, discussed in detail 
bv IDesmet et al. I (|2010l ). The double-peaked H a emission 
line, and the variation of H a and Ha line profiles with the 
orbital phase (lAtwood-Stone et al. IbOlOT ). can be success- 
fully explained with a model of the system with an accretion 
disc around the gainer. IPeters fe Polidan I (J1998D have found 
evidence of disk absorption in the Fe ni ( UV1) resonance 
lines . Moreover , according to the criterion of iLubow fe Shu I 
1975 ), AU Mon should ha ve a permanent accretion disk 



Richards fe Albright Ill999h . 



These spectroscopic indicators and the difficulties in 
reconciling the semidetached model with the observations, 
prompted us to repeat the light curve analysis of C0R0T 
photometry using a binary system model with an accretion 
disk. The model includes active regions on the disk edge, ac- 
counting for the gas stream and the spiral arms in the disk. 
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We aim to demonstrate that the accretion disk model fits 
the CoRoT data better than the semidetached model. 



2 THE MODEL OF THE ACCRETION DISK 

The basic elements of the binary system model with a plane- 
parallel accretion disk and the corresponding light-curve 
synth e sis pr ocedure are described in detail by IDiurasevic I 
(l99lll996h . 

We assume that the disk in AU Mon is optically and 
geometrically thick. The disk edge is approximated by a 
cylindrical surface . In t he current version of the model 
IjDiurasevic et al. I I2OO8I ), the thickness of the disk can 
change linearly with radial distance, allowing the disk to 
take a conical shape (convex, concave or plane-parallel). The 
geometrical properties of the disk are determined by its ra- 
dius (Rd), its thickness at the edge (d e ) and the thickness 
at the center (d c ). This way of approximating the shape 
of the accretion disk is justified by the current hydrody- 
namical mode ling of mass trans fe r in close binary s ystem s 
see, e.g.. iBisikalo et al. N2000I). [Harmanec et al I J2002), 



iNazarenko fe Glazunoval (|2003l , l2006al lbh 

The cylindrical edge of the disk is characterized by its 
temperature, Td, and the conical surface of the disk by a 
radial temperature profile obtain ed by modify ing the tem- 
perature distribution proposed by I Zola I (|199 



% 



T{r) = T d + (T. - T d 



1- 



\Rd-RJ 



(1) 



We assume that the disk is in physical and thermal con- 
tact with the gainer, so the inner radius and temperature of 
the disk are equal to the temperature and radius of the star 
(R,, T»). The temperature of the disk at the edge (Td) and 
the temperature exponent (ot), as well as the radii of the 
star (-R*) and of the disk (-Rd) are free parameters, deter- 
mined by solving the inverse problem. 

The model of the system is refined by introducing active 
regions on the edge of the accretion disk. The active regions 
have higher local temperatures so their inclusion results in 
a non-uniform distribution of radiation. The existence of 
such regions (hot and bright spots) , can be expl ained by the 
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Bisik alo et all (1998. 2000, 2005), Har manec et al 
INazarenko fc Glazunova I (]2003l . l2006al ibTt. Based 



these investigations of systems with accretion disks, we have 
recently updated our model to allow inclusion of up to three 
such active regions: a hot spot (hs) and two bright spots (bsl 
and bs2). Including the active regions in the model leads to 
significant improvement of the quality of the fit. 

The hot spot in our model is a rough approximation of 
the "hot line" which forms at the edge of the gas stream 
flowing from the donor to the disk. Due to the infall of an 
intensive gas-stream, the disk surface in the region of the 
hot spot becomes deformed as the material accumulates at 
the point of impact, producing a local deviation of radiation 
fro m the uniform az imuthal distribution (for more details, 
see IDiurasevic et al. 1120081 ). In the model, this deviation is 
described by the angle # r ad between the line perpendicular to 
the local disk edge surface and the direction of the hot spot's 
maximum radiation. Depending on # ra d, the maximum of 
the hot spot flux can be slightly shifted in orbital phase, 



changing the light-curve asymmetry around the secondary 
maximum and in the region of the primary minimum. 

The bright spots in the model approximate the spiral 
structur e of the disk, pre dicted by hydrodynamical calcu- 
lations (|Heemskerk 1 1 1994T ) . The tidal force exerted by the 
donor star causes a spiral shock, producing one or two ex- 
tended spiral arms in the outer part of the disk. The bright 
spots can also be interpreted as regions where the disk sig- 
nificantly deviates from the circular shape. 



3 THE LIGHT-CURVE ANALYSIS 

This study is based on photometric measurements of AU 
Mon collected by the C0R0T space mission. The data consist 
of 139704 individual measurements with the time sampling 
of 32 seconds, taken during 56 days, encompassing five full 
orbital periods of the system. 

The observ ations were folded to o rbital period using the 
ephemeris from lDesmet et al. I (|2010r ): 

T minl = HJD2454136.6734(2) + ll. d 1130374(l) x E (2) 

Since in this study we were mainly interested in the 
physical and orbital characteristics of the system, we down- 
sampled and divided the data into five light curves, each 
consisting of 1000 points. The light level of the system slowly 
decreases over the time span of the observations, hence we 
normalized each light curve to the maximum of the first one 
at phase 0.25 (the total brightness of the system). In Fig. [T] 
(top five panels), this maximum light level is shown by the 
horizontal dashed line. The normalized light curves were an- 
alyzed separately. 

The light-curve analysis was performed u sing the 
inverse-problem solving method IjDiurasevic lll992T ) based on 
the simplex algorithm, and the above described model of a 
binary system with an accretion disk. 

To obtain reliable estimates of the system parame- 
ters, a good practice is to restrict the number of free pa- 
rameters by fixing them to values obtained from indepen- 
dent sources. Thus we fixed the spectroscopic mass ratio to 
q — M c /Mh — 0.17 (where the subscripts h and c refer to 
the hotter gainer and the cooler donor), and the temper- 
ature of the do nor to T c = 5750Jf , based on the spectro- 
scopic study bv lDesmet et al. I (|2010n . The temperature of 
the gainer was a free parameter in the preliminary analysis, 
and for the final solution we adopted the average of the indi- 
vidual solutions for the five successive C0R0T light curves. 
The average value of Th — 15870A', obtained in this way, 
was then kept fixed in the final analysis. We adopted the 
factor of non-sync hronous rotation of the gainer //, = 5.2 
from the study bv ldazunova et al. I (|2008l ). and considered 
the rotation of the donor as synchronous (/ c = 1.0), assum- 
ing that it has filled its Roche lobe (i. e. the filling factor 
of the donor was set to F c = 1.0). In addition, we set the 
gravity darkening coefficient and the albedo of the gainer 
to Ph — 0.25 and Ah = 1.0 in accordance with von Zeipel's 
law for radiative shells and complete re-radiation; for the 
donor we set /3 C — 0.08 and A c = 0.5, as is a ppropr i ate fo r 
stars with conve ctive envelopes acc ording to iLucv I |l967J), 
iRucinski I j 19691 ) and lRafertl |l980T ). 

The limb-darkening for the components was taken into 
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accou nt by using the non-linear approximation by Idaretl 
(J200CJ), which can be expressed as 



^ = l-a 1 (l- M 1 / 2 )-a 2 (l- M )-a 3 (l-M 3/2 )-a4(l-^ 2 ) 

where 01,2,3,4 are the limb-darkening passband-specific co- 
efficients, 1(1) is the passband-specific intensity at the cen- 
ter of the stellar disc, and fi = cos 7, where 7 is the an- 
gle between the line of sight and the local surface normal. 
The central intensity is calculated for the effective wave- 
lengths of the observations (600 nm for the C0R0T light 
curves and 550 nm for the V-band light curve), using a sim- 
ple blackbody approximation. For a given metallicity, the 
values of the passband-specific limb-darkening coefficients 
are derived from the current values of the stellar effective 
temperature T e fr and surf ace gravity log g in each iteration, 
by bi-hnear interpo lation (Press et al. 1 11992). of both quan- 
tit ies from tab les of lClaretl ( 20000 for the V light curves and 
of ISing I l|2010h for the C0R0 T light curves. The proc edure 
is described in more detail bv lDiurasevic et al. I 1)20041 ) . The 
limb-darkening was applied to the disk in the same way, with 
log g corresponding to the middle of the disk radius. 



4 RESULTS 

The results of the light-curve analysis based on the described 
model of AU Mon are given in Table [T] The first column con- 
tains parameter designations, and the following five columns 
give the values derived from each of the five C0R0T light- 
curves, with the mean values and their estimated uncertain- 
ties. The uncertainties were estimated from a set of solutions 
obtained for the five observed light curves and three different 
values of the mass ratio: q\ = 0.14, qi = 0.17 and 53 = 0.20 
(c hosen according to t he error assigned to the mass ratio 
bv lDesmet et al. II2OIOI . q = 0.17 ± 0.03), resulting in a to- 
tal of 15 values for each parameter. The uncertainties given 
in Table [l] are the maximal deviations of these values from 
the mean. Table [J als o lists the results of applying a simple 
Roche model (see e.g. iDiurasevic | [l992) to the C0R0T light 
curves and the results of applying the accretion disk model 
to the ground-based V-band light curves (discussed in de- 
tail in Section [4. ip . The first three rows of Table [T] present 
the number of points in the light curve (n), the final sum 
of the squares of the residuals between the observed (LCO) 
and the synthetic (LCF) light curves, J^(0 — C) 2 , and the 
root-mean-square of the residuals <x rms , respectively. 

The best fit model of AU Mon contains an optically and 
geometrically thick accretion disc around the hotter, more 
massive gainer star. With a radius of Rd w 13Rq , the disk is 
more than twice as large as the central star (Rh ~ 5Rq). The 
disk has a moderately concave shape, with central thickness 
of d c ~ OARq and the thickness at the edge of d e fa 1.6Rq. 
The temperature of the disk increases from Td = 51907^ at 
its edge, to Tu = 15870-ft' at the inner radius (where it is in 
thermal and physical contact with the gainer), according to 
Eq. [l] with the temperature profile exponent ar — 6.5. The 
effective temperature of the disk is significantly higher than 
the temperature at its edge. 

We were able to model the asymmetry of the light curve 
very precisely by incorporating three regions of enhanced 



radiation on the accretion disk: the hot spot (hs), and two 
bright spots (bsl and bs2). The hot spot (hs) is situated 
at longitude Xhs ~ 330°, roughly between the components 
of the system, at the place where the gas stream falls onto 
the disk. The longitude A is measured clockwise (as viewed 
from the direction of the +Z-axis, which is orthogonal to 
the orbital plane) with respect to the line connecting the 
star centers (+X-axis), in the range 0° — 360°. The tem- 
perature of the hot spot is approximately 70% higher then 
the disk edge temperature, i.e. Ths ~ &700K. The hot spot 
can be interpreted as a rough approximation of the "hot 
line" which forms at the edg e of the gas stream betwee n 
the components. According to lAtwood-Stone et al. I (|2010T ) . 
a gas-stream with a temperature near 8000A' can explain 
the excess red- and blue-shifted H a emission near phases 
0.2 and 0.7, respectively. We note that this spectroscopic re- 
sult is an independent confirmation of our photometrically 
estimated temperature of the gas stream. 

Although including the hot spot region into the model 
significantly improves the fit, it cannot explain the light- 
curve asymmetry completely. By introducing two additional 
bright spots (bsl and bs2), larger than the hot spot and 
located on the disk edge at \b a \ ~ 170° and A& S 2 fa 50°, the 
fit becomes much better. 

The bright spots can be related to the spiral shocks that 
result from radiative cooling and form at the outer boundary 
of the disk. Since the disk is large, filling about 56% of the 
gainer's critical Roche surface, the tidal forces exerted by 
the don or can cause a spir al-shaped tidal shock in the disk - 
see e.g. iHeemskerk I |l994J). Such a shock wave can produce 
one or two extended spiral arms in the outer parts of the 
disk. 

The first arm, represented in our model by a bright spot 
(bsl), is located on the disk edge at longitude Ab s i fa 170°. 
This is also a region where we can expect loss of matter from 
the gas stream and the disk through the Lagrangian point 
L3, forming some kind of a circumbinary shell. 

The fit was additionally improved by introducing the 
second bright spot (bs2), located at longitude Xbs2 ~ 50°. 
This spot is the largest active region, with a temperature 
about 30% higher than the disk edge temperature. It can be 
interpreted as the second spiral arm in the disk. 

We note that the system can also be modeled with ac- 
tive regions (dark spots) on the donor, and without the ac- 
tive regions on the accretion disk. Such a model would ex- 
plain the period-to-period variations in the light curves by 
the presence, development and migration of spots over the 
surface of the donor. However, the model with active re- 
gions on the accretion disk seems to be more appropriate. 
Namely, the relatively fast variations of the light curves are 
more likely to originate from the changes in the disk struc- 
ture, produced by variable mass outflow from the donor, 
than from the motion of stellar spots, not expected to take 
place on these timescales. 

In order to make a comparison betwee n the results of 
our study and that of iDesmet et al. I 1)20101 ). we made sev- 
eral trial runs with a simple semidetached Roche model of 
AU Mon. The semidetached model cannot fit the obser- 
vations as well as the model with an accretion disk. The 
fit can be improved by using an anomalousl y large gravity 
darke ning exponent of the donor, as done bv lDesmet et al. I 
1)20101 ). Our previous investigations of anomalously high val- 
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Table 1. Results of the analysis of CoRoT and V-band light-curves, obtained by solving the inverse problem for the Roche model with 
an accretion disk around the more-massive (hotter) component, and with the simple semidetached Roche model 



Quantity 


CoRoT (1) 


CoRoT (2) 


CoRoT (3) 


CoRoT (4) 


CoRoT (5) 


MEAN CoRoT 


MEANcoRoT 


V - filter 




Disk model 


Disk model 


Disk model 


Disk model 


Disk model 


Disk model 


Roche model 


Disk model 


n 


1000 


1000 


1000 


1000 


1000 






2682 


E(O-C) 2 


0.0497 


0.0410 


0.0636 


0.0480 


0.0454 






2.6880 


^rms 


0.0070 


0.0064 


0.0080 


0.0069 


0.0067 






0.0317 


q 


0.17 


0.17 


0.17 


0.17 


0.17 


0.17 ±0.03 


0.17 ±0.03 


0.17 


i[°] 


80.07 


80.09 


80.06 


80.09 


80.07 


80.1 ±0.6 


78.8 ± 0.6 


80.08 


F d 


0.575 


0.578 


0.532 


0.560 


0.518 


0.55 ±0.03 




0.575 


T d [K] 


5140 


5174 


5209 


5256 


5192 


5190 ± 70 




5250 


de [a orb ] 


0.040 


0.039 


0.039 


0.039 


0.039 


0.040 ±0.005 




0.042 


dc [a orb ] 


0.010 


0.010 


0.010 


0.010 


0.010 


0.01 




0.01 


a T 


6.50 


6.50 


6.50 


6.50 


6.50 


6.5 




6.5 


F h 


0.565 


0.565 


0.565 


0.565 


0.565 


0.565 ± 0.002 


0.62 ±0.01 


0.565 


Th[K] 


15894 


15905 


15886 


15902 


15890 


15890 ± 400 


14100 ± 400 


15890 


Tc[K] 


5750 


5750 


5750 


5750 


5750 


5750 


5750 


5750 


Ahs = Thg/Tj 


1.68 


1.68 


1.70 


1.70 


1.68 


1.7 ±0.1 




1.68 


hs [°] 


19.8 


19.5 


19.7 


19.9 


19.6 


19.7 ±0.4 




20.0 


A hs [°] 


325.6 


329.0 


330.7 


334.2 


337.6 


331 ±6 




312.4 


^rad[°] 


-33.2 


-34.8 


-34.6 


-34.7 


-32.1 


-34 ±7 




-32.1 


Absl = T bsl /T d 


1.34 


1.47 


1.35 


1.48 


1.51 


1.4 ±0.2 




1.31 


O h Bl[°] 


39.4 


39.6 


39.3 


30.7 


17.0 


33 ±12 




36.5 


^bsl[°] 


177.1 


178.5 


178.4 


164.5 


158.0 


171 ± 13 




117.9 


Abs2 = T bs2 /Td 


1.39 


1.41 


1.28 


1.33 


1.25 


1.3 ±0.1 




1.50 


^bs2[°] 


49.7 


49.4 


49.2 


49.8 


48.5 


49 ±6 




50.0 


^bs2 [°] 


60.0 


47.2 


52.6 


51.0 


53.8 


53 ±14 




73.7 


n h 


8.614 


8.613 


8.611 


8.614 


8.616 


8.6 ±0.5 


7.9 ±0.6 


8.612 


H c 


2.156 


2.156 


2.156 


2.156 


2.156 


2.2 ±0.2 


2.16 ±0.05 


2.156 


M h [M@] 


6.957 


6.957 


6.957 


6.957 


6.957 


7.0 ±0.3 


7.0 ±0.4 


6.957 


M C [M@] 


1.183 


1.182 


1.183 


1.182 


1.182 


1.20 ±0.2 


1.2 ±0.2 


1.183 


^hPo] 


5.081 


5.082 


5.083 


5.081 


5.080 


5.1 ±0.5 


5.6 ±0.5 


5.082 


^ C [R ] 


10.083 


10.083 


10.083 


10.083 


10.083 


10.1 ±0.5 


10.1 ±0.5 


10.083 


log<? h 


3.868 


3.868 


3.868 


3.868 


3.869 


3.87 ±0.05 


3.78 ±0.06 


3.868 


log 9c 


2.504 


2.504 


2.504 


2.503 


2.503 


2.50 ±0.02 


2.50 ±0.02 


2.504 


bol 


-3.138 


-3.141 


-3.136 


-3.140 


-3.136 


-3.1 ±0.4 


-2.8 ±0.4 


-3.137 


M bol 


-0.210 


-0.210 


-0.210 


-0.210 


-0.210 


-0.21 ±0.09 


-0.21 ±0.09 


-0.210 


a orb [R©] 


42.11 


42.11 


42.11 


42.10 


42.10 


42.1 ±0.4 


42.1 ±0.5 


42.1 


n d [Ko\ 


13.20 


13.28 


12.23 


12.87 


11.90 


12.7 ±0.6 




13.2 


de[R ] 


1.69 


1.66 


1.65 


1.66 


1.66 


1.66 ±0.05 




1.78 


dc[R Q ] 


0.42 


0.42 


0.42 


0.42 


0.42 


0.42 ±0.01 




0.42 



FIXED PARAMETERS: q = A4 c /Mh = 0.17 - mass ratio of the components, T b = 15890K - temperature of the more-massive (hotter) 
gainer, T c = 5750X - temperature of the less-massive (cooler) donor, F c = 1.0 - filling factor for the critical Roche lobe of the donor, 
/ b = 5.2; f c = 1.00 - non-synchronous rotation coefficients of the gainer and donor respectively, /3 h = 0.25; /3 C = 0.08 - gravity-darkening 
coefficients of the gainer and donor, A b = 1.0; A c = 0.5 - albedo coefficients of the gainer and donor. 

QUANTITIES: n - number of observations, S(0 — C) 2 - final sum of squares of residuals between observed (LCO) and synthetic (LCF) 
light-curves, cr rms - root-mean-squarc of the residuals, i - orbit inclination (in arc degrees), Fd = Rd/Ryc - disk dimension factor (ratio 
of the disk radius to the critical Roche lobe radius along y-axis), Td - disk-edge temperature, d c , d c , - disk thicknesses (at the edge 
and at the center of the disk, respectively) in the units of the distance between the components, ax - disk temperature distribution 
coefficient, F b = Rh/Rzc - filling factor for the critical Roche lobe of the hotter, more-massive gainer (ratio of the stellar polar radius 
to the critical Roche lobe radius along z-axis), A bs bsl.bs2 = T bs bsl,bs2/Td - hot and bright spots' temperature coefficients, # bs bs i bs 2 
and A bs bs i bs 2 - spots' angular dimensions and longitudes (in arc degrees), <9 ra d - angle between the line perpendicular to the local 
disk edge surface and the direction of the hot-spot maximum radiation, Q b c - dimensionless surface potentials of the hotter and cooler 
components, M h c [Mq], 1Z h C [R0] - stellar masses and mean radii of stars in solar units, log g h - logarithm (base 10) of the system 



components effective gravity, M' c . - absolute stellar bolometric magnitudes, a orb [Rq], KjJRq], d e [RQ], d c [Ro] - orbital semi-major 
axis, disk radius and disk thicknesses at its edge and center, respectively, given in the solar radius units. 
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Figure 1. Left: CoRoT observed (LCO) and synthetic (LCC) light-curves of AU Mon (the normalized maximum light level is shown by 
the horizontal dashed line) obtained by applying a simple semidetached Roche model (dashed line) and an accretion disk model (solid 
line) with the active regions on the disk edge. The successive individual light curves are shifted by 0.2 magnitude; Right: the final O-C 
residuals between the CoRoT observed (LCO) and optimal synthetic (LCC) light curve, obtained by applying the simple semidetached 
Roche model (Roche) and accretion disk model (Disk). 
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Figure 2. The views of the optimal accretion disk model of the system AU Mon at orbital phases 0.00, 0.30, 0.47 and 0.70, obtained 
with parameters estimated by the CoRoT light curve analysis. 



6 G. Djurasevic et al. 



1.0 


X 




DONOR 


AU 


Mon 




CoRoT 


0.8 


3 














0.6 


Z 








2 

3 






0.4 
0.2 


J:*'"'" 


"""V 


DISK^^. 




^^i* 


^iV 


fcZv^^Z 


0.0 
0. 














^ PHASE 


0.1 


0.2 


0.3 0.4 


0.5 0.6 


0.7 O.B 


0.9 


1.0 1.1 1.2 1.3 



Figure 3. The calculated fluxes of the donor and of the accretion 
disk, normalized to the donor flux at the phase 0.25 and given for 
all five CoRoT light curves. 



ues of gravity darkening exp onents reported in literature 
IjDiurasevic et al. 1120 03. 2006) have shown that such anoma- 
lies are usually the consequence of an inadequate model used 
for the light curve analysis. The parameters we obtained for 
a semidetached Roche model (without the disk), with the 
theoretical value of gravity darkening are shown in Table [l] 
(column 7) . The semidetached model gives a slightly smaller 
inclination (by approximately 2%), and a larger radius of 
the gainer (by approximately 9%), with a lower temperature 
(by approximately 10%), which was to be expected since the 
large accretion disk, surrounding the gainer, is characterized 
by a lower effective temperature. In summary, the absence 
of the accretion disk in the semidetached model is compen- 
sated by a smaller inclination, and a larger, cooler gainer. 
Since the s pectral type of the gainer was reported to be 
B5 (see e.g. lDesmet et al. II20101 . and the references within), 
the gainer's temperature estimated with the accretion disk 
model is in better agreement with this spectral type than the 
temperature estimated with the simple semidetached model. 

Another prob le m wit h the semidetached model used 
by iDesmet et al. I (|2010r ) is its inabilit y to rep r oduce 
the Rossiter-McL aughlin (RM) effect iJRossiter I 1 1924 
iMcLaughlin Ill924f ). To improve the q uality of the ra d ial ve - 
locity curve fit to the observations, IDesmet et al. I (|2010r ) 
had to enlarge the radius of the donor much over the value 
obtained from the light curve analysis (see Figure 11 of their 
paper). This discrepancy can also be solved by introduc- 
ing an accretion disk around the gainer. Namely, from the 
bottom panel of Fig. [4] which shows the system at orbital 
phase 0.5 (in the secondary eclipse), it is intuitively clear 
that the accretion disk contributes to the RM effect in a 
similar way as the enlargement of the donor. This remains 
to be proven by spectroscopic modeling. Nevertheless, our 
qualitative conclusion is that the enlargement of the donor 
is unnecessary if the system is modeled with an accretion 
disk. 

The fits are shown in the left panel of Fig. Q] for each 
orbital period. The successive individual light curves are 
shifted by 0.2 mag to save space. The figure shows the obser- 
vations (LCO - CoRoT) and two fits: the fit of the described 
accretion disk model (LCC - Disk) and the fit of a simple 
semidetached Roche model (LCC - Roche) with the fixed 
gravity darkening exponent. The residuals for both models 
are given in the right panel of Fig. [T] The model with the 
accretion disk, made with the theoretically expected grav- 



ity darkening exponent j3 = 0.08 for the convective donor, 
evidently provides a superior fit to the observations. 

Figure [2] shows the appearance of AU Mon in orbital 
phases 0.00, 0.30, 0.47 and 0.70. The phases were chosen so 
as to make both components and the active regions on the 
disk edge visible. 

Finally, Fig. [3] presents the individual synthetic fluxes 
of the donor and the disk, normalized to donor brightness at 
orbital phase 0.25. The main contribution to the flux comes 
from the hotter, more massive gainer, and is approximately 
3 times greater than the contribution from the donor, so it 
is not shown here to save space. In the present model, the 
donor and the gainer fluxes are the same for all five orbital 
periods in CoRoT data, and the shape variations of the light 
curves are entirely due to the changes in the flux emitted by 
the accretion disk, probably produced by the variable mass 
outflow from the donor. 



4.1 The long-term variability 

The long-term variability in total brightness of the system 
was firs t explai n ed as the re sult of fluctuating mass transfer 
rate bv lPeters I (|l99ll . ll994T ). The rate of mass transfer may 
change periodically due to the oscillations of the donor star 
around its critical Roche lobe, causing changes in the size 
and density of the accretion disk and/or circumbinary mat- 
ter. IVivekananda fc Sarma | (119981) . who s tudied the pho- 



toelectric measurements of 



Lorenzi (Il980bf) . using a similar 



semidetached model as lDesmet et al. 1 ( 20101 ). suggested that 



the long-term variation may be due to the precession of the 
accretion disk around the gainer, caused by the gravitational 
perturbations by the donor. 

The ephemer i s of t he long-term variability calculated 
by IDesmet et al. I (|2010l ) is: 



= HJD2443105.1(±1.4) +416. 9(±8.7) x E. 



(3) 



This puts the CoRoT observations at the minimum of 
the total light, between phases 0.46 and 0.59 of the long-term 
variation. 

In order to study the behavior of the model with accre- 
tion disk during the entire lon g cycle, we used the gr ound 
based photometric da ta from lLorenzi I l|l980bl . 119851 ) and 
iKilkennv et al. I (|1985l ) in the V passband (courtesy of M. 
Desmet). 

After folding the light curves to the long period accord- 
ing to the ephemeris given by Eq. [3] we attempted to model 
the data from the maximum (phases 0.9 to 1.1) and from 
the minimum (phases 0.4 to 0.6) of the long period, using 
the model obtained with CoRoT data, with the quantities 
related to the accretion disk as free parameters. We found 
that the variation in total brightness cannot be explained 
by changes in the size of the accretion disk and the sizes 
and locations of the active regions on the disk edge. Such 
changes would cause variations in depth of the light curve 
minima, but there is no evidence that the amplitude of the 
light curve is changing over th e long period. Our an alysis 
thus confirms the conclusion of IDesmet et al. I (J2010r) . that 
the long-term variation must be the result of variable atten- 
uation of the total light of the system by some form of a 
circumbinary shell. 

The accretion disk model, when used with the same 
methods and free parameters we employed in the analysis of 
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Figure 4. From top to bottom: V-band observed light curves 
(LCO) at maximum (max) and minimum (min) of the long-term 
period (Pi ong ~ 417 days) and corresponding synthetic (LCF) 
light-curves of AU Mon obtained by applying accretion disk model 
with active regions (AD+AR) - the normalized maximum light 
level is shown by the dashed line; The final O-C residuals for the 
maximum and minimum light curves obtained by applying accre- 
tion disk model with active regions (AD+AR); Observed (LCO) 
V-band light curve (after removal of the long-term period) and 
the corresponding fitting curves (LCFad+ARi LCFad) obtained 
by applying accretion disk model with (AD+AR) and without 
(AD) active regions; the corresponding O-C residuals obtained 
for the accretion disk model with active regions (AD+AR); the 
view of the model of the system AU Mon at orbital phase 0.50 
obtained with parameters estimated from the V-band light curve 
analysis. 



the CoRoT data, fits the ground-based data very well and 
the parameters obtained in this way (given in column 8 of 
TableLT} are in good agreement with those obtained from the 
CoRoT light curves. We conclude that our model describes 
both data sets adequately. 

The ground-based data and the synthetic light-curves 
from our model are shown in Figure [4] The top panel shows 
the fit to the light curves in the minimum and maximum of 
the long cycle, the middle panel shows the fit to the entire 
ground-based dataset, reduced for the effects of the long- 
term variation, and the bottom panel gives a representation 
of the model obtained from the entire ground-based dataset 
in the phase of the secondary eclipse. 

4.2 The short-term variability 

There are evident short-term variations in the CoRoT light 
curves of AU Mon. The residuals after subtracting the syn- 
thetic light curve from the observations (Figure [T] right 
panel), also show quasi-periodic variability. Two features 
stand out: a lump in the secondary minimum (around phase 
0.5) which is present to some degree in all five residual curves 
and is partially due to the inability of the model to fit the 
observations perfectly. However, the shape and amplitude of 
the lump change from period to period, suggesting the pres- 
ence of intrinsic variability in either the G type secondary 
or in the accretion disk or both. 

The other notable feature in the residual light curves is 
the fast and seemingly periodic variation inside the p rimary 
minimum (around phase 1.0). iDesmet et al. I ([2010), who 
also noted the fast variation, suggested that its origin is a 
non-uniform brightness distribution of the accretion disk. 

There is, however, an alternative interpretation: we be- 
lieve that the fast variation during the primary minimum of 
AU Mon is caus ed by a hidden mode o f oscillation of the 
primary (see e. g. iMkrtichian et al. 1120051 ). The origin of the 
variation is unlikely to be the secondary star or the accretion 
disk, because these components are at least partially visible 
during the entire orbital cycle, while the fast variation oc- 
curs only in the primary minimum. If our interpretation is 
correct, it also represents an indication that the primary star 
is indeed producing oscillations. 

We performed a frequency analysis of the residuals re- 
maining after fitting our model to the observations using Pe- 
RIOD04, a sign al analysis program b ased on Fourier Trans- 
form methods (|Lenz fc Breger II2005T ). A time series of 5000 
points was formed by concatenating the residuals remaining 
after fitting each orbital period, allowing the computation 
of a periodogram in the frequency range from to 45 gT 1 
with a frequency resolution of 0.0001 dT 1 , shown in part in 
the upper panel of Figure LU 

The low-frequency ra nge is dominat e d by the harmon- 
ics of the orbital period. IDesmet et al. I (|2010l ) were faced 
with the same problem during the frequency analysis of their 
residual data, and dealt with it first by cutting out the pri- 
mary minima from their data set, and then by dividing the 
data with a spline-fit through the highest peaks in the PDM 
periodogram, with an effect of a high-pass filter. 

Our attempt to remove the harmonics of the orbital 
frequency by cutting out regions around the primary and 
secondary minima did not result in a significant improve- 
ment in the appearance of the periodogram, so we used 
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Figure 5. The Fourier periodogram computed from the O — C 
residuals in the frequency range from to 12 d , before (up- 
per panel) and after (lower panel) prewhitening the data for the 
harmonics of the orbital frequency. 



Table 2. Frequencies found in the O — C residuals. 

ID Frequency Amplitude SNR 

d~ 1 fiHz mmag 

/l 10.4079 120.46 0.28 7.5 

h 8.3036 96.11 0.21 5.3 



the method of iterative sine fitting (prewhitening) to re- 
move about 40 harmonics from the signal. The resulting, 
cleaned periodogram is shown in the lower panel of Figure [5] 
The arrows indicate the two stable frequencies reported by 
iDesmet et al~l l|2010l ). at 10.40 and 8.30 d' 1 . The frequen- 
cies and amplitudes of these two peaks, as determined from 
our residual data (listed in Table [2]), are in good agreement 
with their results, indicating that in this case the frequency 
analysis is not very sensitive to the model of the binary sys- 
tem. 



5 SUMMARY 

The motivation for this work came from the recent spectro- 
scopic studies that offer compelling evidence of the presence 
of an accretion disk in the interacting eclipsing binary AU 
Mon. We analyzed the CoRoT photometric observations of 
AU Mon using a model that includes a geometrically and 
optically thick accretion disk. The model fits the observa- 
tions very well. We note that we also modeled the system 



using a simple Roche model without an accretion disk, and 
found the final sum of the squares of the residuals between 
the observed and synthetic light-curves for this model to be 
significantly larger than for the model with the disk. 

The accretion disk in our model is moderately con- 
cave and extends to roughly 60% of the critical Roche lobe 
of the gainer. The asymmetry of the light curve is suc- 
cessfully modeled by incorporating three active regions on 
the accretion disk: a hot spot (hs), and two bright spots 
(bsl and bs2), which in turn represent the "hot line" and 
two spiral arms of the disk and/or the deviations of the 
disk from a circular shape. The temperature of the hot 
spot (Ths ~ 8700if) is in good agreement with the gas- 
stream temperature T ~ 8000i(" obtained from H a emis- 
sion line profile modeling n ear orbital phases 0.2 and 0.7 
l|Atwood-Stone et al. Il2010h . which is an independent con- 
firmation of our ph otometrically estimat e d tem perature of 
the hot line region. lAtwood-Stone et al. I (J2010j) found that 
a disk with a maximum temperature of T ma x ~ 13700-K" is 
in general agreement with the observed UV spectrum. The 
estimated temperature of our disk is in the range from 5100 
K (at the edge of the disk) to 15870 K (in the inner part 
that is in contact with the gainer). The radial temperature 
profile of the disk, given by Eq. [T] sho ws that the maximum 
tempe rature of the disk estimated by lAtwood-Stone et al. I 
(J2010j) lies well within the range covered by our model. 

The period-to-period variability of the CoRoT light- 
curves can be explained by the changes in the contribution 
of the disk flux to the total light, due to variations in the 
outflow from the donor. However, the long-term variability 
cannot be explained by the changes in the disk parameters 
alone, because such changes would cause variations in the 
amplitude of the light curve over the long period. As the 
changes of the light curve amplitude were not observed, the 
long-term variability of the system can only be interpreted 
as a result of a periodic change in the attenuation of the 
total light by an inhomogeneous circumbinary shell. 
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